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ABSTRACT 


A direction finding antenna based on residue number systems (RNSs) is presented. 
Special spacing requirements for the array elements are derived and processing of the array 
output for high speed direction of arrival (DO A) estimates is discussed. 

The RNS antenna processor encodes the interferometer phase response, which is a 
sawtooth folding waveform. By design, the phase response of each element pair folds with 
folding period equal to the chosen modulus. 

If DOA samples are generated by stepping the emitter direction between -90 and 90 
degree in small increments, some samples will fall about the code transition points and 
result in encoding errors. The resulting DOA estimates from such a realizable system could 
contain large spikes or “glitches” at these points. These encoding errors in the resolved 
DOA can be reduced by interpolation. Three primary methods will be discussed to compare 
the capability of removing the glitches: LSB-Shift Method, Random-LSB Method, and 
Shift Last Good-Sample Method. A comparison of the performance of the three methods is 
made on the basis of simulation data. 
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I. 


INTRODUCTION 


A. OVERVIEW 

Direction fin din g (DF) is the area of electronic-warfare support (ES) whose 
objective is to find the direction of arrival (DOA) of a targeted transmitter. ESM performs 
search, interception, location, and identification of hostile radiations. DF systems passively 
examine the spectrum use by hostile emitters and process the signals to obtain enemy DOA 
information. 

A general assumption in most DF systems is that the received field exhibits far-field 
plane-wave behavior with linear polarization. In reality, incident fields are often non-planar 
with phase-front distortion caused by multipath and scattering. Typical DF system 
techniques assume that only a single signal source is received. Because of multipath 
ionospheric propagation (skywave), and wave scattering by objects close to the DF site or 
co-channel transmitting stations, signals via multiple paths and sources are often received at 
DF sites. A DF system designed to handle one source at a time can work in a multiple 
source environment if other signals are well separated in frequency, time, or direction of 
arrival. 

In practice, DF stations often receive signals from several sources in the same 
frequency range at the same time, and from many different angles of arrival [Ref. 1]. To 
resolve multicomponent wavefields, modem techniques of spectral estimation methods are 
used. These methods include eigenanalysis and linear prediction. 

Modem direction finding systems fall into three main categories: interferometer- 
techniques, scanning beam, and simultaneous-multiple beam. Interferometer systems based 
on phase-comparison techniques, have the advantage of a fast response, but generally use 
wide coverage antennas which result in low sensitivity. In addition, they require relatively 
complex microwave circuitry that must maintain a precise phase match over a wide 
frequency band under extreme environmental conditions. When high accuracy is required, 
wide baseline interferometers are used with ambiguity resolving circuitry. The mechanically 
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sc annin g beam requires only a single antenna. However, this technique has the disadvantage 
of a low probability of intercept. The simultaneous multiple-beam system uses an antenna 
or several antennas, forming a number of simultaneous beams, thereby retaining the high 
sensitivity of the scanning beam approach while providing fast response. However, it 
requires many parallel receiving channels, each with full frequency coverage. 

This thesis is concerned with a technique for radiating element placement and signal 
processing based on the residue number systems (RNS). In principle, the direction of arrival 
can vary in azimu th (0-360 degrees) and elevation (0-90 degrees). However, this paper 
deals exclusively with the azimuth angle of arrival for a linear array. In addition, an 
underlying assumption is that the receiver is far enough from the transmitter so that the 
arriving wavefront is essentially planar. 

The basic concept of residue arrays' was presented in a previous thesis by Rodrigues 
[Ref. 2]. A simple block diagram of such a system is shown in Figure 1. It is essentially an 
interferometer type DF system consisting of N -1 pairs of elements. Each pair includes 
element 1 so that there are N -1 pairs if the total number of elements is N. The number of 
pairs is determined by the number of moduli employed in processing the signal. It will be 
seen that the number of moduli controls the angle resolution, the element spacings, and the 
number of comparators required in the processor. The phase difference between each pair 
of elements is obtained using a phase detector. In principle, phase or amplitude comparison 
can be used, but Rodrigues has shown that the phase is less sensitive to variations in the 
DOA and element factors. The phase differences are quantized and encoded using the RNS 
[Ref. 3]. The encoded data from all pairs is combined to obtain a DOA estimate. 

This DF method has several attractive features. First, arbitrarily fine unambiguous 
angle resolution can be achieved with as few as three elements. Secondly, the method is 
fast, providing a DOA estimate limited only by the RNS processor speed. The processor 


1 This term will be used throughout the thesis to refer to arrays with element locations or 
encoding algorithms based on residue or symmetrical number systems. 
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employed here is based on parallel digitization of the N -1 phase signals, and therefore 
nearly ins tantan eous DOA estimates are obtained. 



DOA 


Figure 1. DF Array with RNS Processing 

The RNS DF method has several potential problems. They are primarily 
performance limitations imposed by imperfections in the system hardware. Several of these 
occur in the array beamforming network: 

1. Mutual coupling variations due to the non-periodic separation of elements. 

2. Imbalance in the weights and transmission lines. 
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3. Low accuracy of the phase comparators due to low signal levels. 

Other errors are introduced in the RNS processing: 

1. Encoding errors that arise from variations in the comparator thresholds. 

2. Increase in quantization error due to frequency deviations. 

The most serious problem is the encoding errors which, because of the non-sequential 
binary encoder outputs, can cause extremely large errors in the DOA estimate. In fact, the 
net effect of all of the errors listed above can be reduced to an equivalent encoding error. 

This thesis examines three interpolation schemes that can be used to reduce the 
encoding errors which result from quantization in the RNS signal processing. Simulation 
d ata shows, for example, that small errors in the comparator thresholds can result in large 
DOA errors. Interpolation methods can be used to reduce the occurrence and magnitude of 
the errors. 

Three types of interpolation are investigated: (1) a shifted least significant bit (LSB) 
method, (2) a randomized LSB method, and (3) a shift last good sample method. The 
effectiveness of these three approaches is examined using computer simulation. 

B. THESIS OUTLINE 

Chapter II summarizes common direction finding techniques such as amplitude 
response, phase difference, and time delay methods. Chapter III covers array beamforming 
and the mathematical development for a linear array antenna. Also, direction finding for 
multiple emitter DOAs is introduced by using the array correlation matrix. Chapter IV is 
devoted to the process of encoding the phase response with the RNS. First the RNS 
encoding procedure is briefly introduced. The RNS antenna architecture is then derived 
with simulation results shown to demonstrate the array performance. Next, we examine 
encoding the phase response with the RNS. The RNS encoding procedure is briefly 
described and simulation results are presented. In Chapter V, the encoding error correction 
algorithms are derived. Three primary methods are investigated to remove the quantizing 
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errors or “glitches.” A comparison of the performance of the three methods is made on the 
basis of simulation data. Chapter VI presents a summary, some concluding remarks, and 
recommendations for future work. 
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H. DIRECTION FINDING PRINCIPLES AND METHODS 


As early as 1899 it was recognized that electromagnetic waves could be aimed or 
beamed in a specific direction by the use of antennas such as loops or spaced dipoles. For 
receiving purposes this meant that the direction of arrival of a radio wave could be 
determined by moving or steering a directive antenna beam in azimuth until a maximum 
signal was obtained. The beam could be steered by physically moving an antenna or by 
appropriately summing the outputs of the individual elements of an array. 

The science of direction finding has been the subject of much study since the early 
1900’s. A major improvement over early DF systems came with the development of the 
Adcock DF system employing four or more spaced elements with an overall aperture of 
one-half wavelength or less. The Adcock system achieves an instrumental accuracy of about 
5 degrees. 

Since the end of World War II, there has been a trend toward larger aperture 
systems. For example, Doppler DF systems, in which frequency modulation (FM) by a 
sinusoid is induced on the received signal by a rotating antenna or its equivalent, usually 
range in size from one-half to five wavelengths and are capable of accuracies of 2 degrees 
or less [Ref. 4], 

The Wullenweber DF system, which was developed in Germany in the early 1940's, 
consists of a large circular array on the order of ten wavelengths in diameter [Ref. 5]. At any 
instant of time only one segment of the array is used, and with proper phasing of the 
elements in use, this segment is made equivalent to a planar array. By commutation the 
planar array is in essence caused to move in a circle. An accuracy of 0.5 degree or less is 
possible with this system. 

At mi crowave frequencies, antennas having aperture sizes of several wavelengths 
and with extremely directive patterns can be physically rotated to give accurate DOA 
information. An alternate method which is very accurate for large values of signal-to-noise 
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ratio (SNR) is the amplitude-sensing monopulse technique in which two identical directive 
antennas are placed some distance apart. 

A. SPATIAL DIRECTION-FINDING PRINCIPLES 

The main function of a direction finder is to find out the direction of arrival of an 
incident electromagnetic wave relative to the coordinates of the direction finding site. 
Figure 2 shows a representative direction-finding spatial coordinate system with the DF 
located at the origin coordinate point. The direction of arrival is specified by the azimuth 
angle, measured from the y-axis, and the elevation angle measured from the z-axis. 

General direction- finding principles assume that the electromagnetic field incident 
on the DF exhibits a far-field planar wave structure with linear polarization. As depicted in 
Figure 2, the radial E-field is zero. The E-field, and the H-field are in space quadrature 
[Ref. 6], 



The direction of propagation is indicated by Poynting's vector, P . In practice, the incident 
electromagnetic field is usually nonplanar with phase-front distortion, created by multipath 
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and scattering, and depolarization produced by nonuniform ionospheric propagation 
effects. A full-capacity generic DF should measure direction of arrival in three-dimensional 
space. However, many operational situations require only measurement of the directional 
component in the azimuth plane. 

The measured angle of arrival, <(i^, is usually called a bearing angle or line of 
bearing (LOB). Before examining residue arrays, the capabilities of several conventional 
DF methods are presented. 

B. DIRECTION FINDING TECHNIQUES 

1. Overview 

Common direction-finders determine direction-of-arrival information by one of 
three measurement methods [Ref. 6]: 

1. Amplitude response, 

2. Phase difference, 

3. Time delay. 

The angle of arrival is converted into a voltage analog using amplitude response, phase 
difference, or time delay. Basic conversion techniques are presented in the following 
sections. 

2 . Amplitude Response 

Two types of amplitude response methods are used to obtain DF information: direct 
and comparative. For direct amplitude response, the directional properties of elemental H- 
field loops and E-field dipoles provide response minima as they are rotated in azimuth. 
Bearing data are acquired at these minima response orientations. 

Elemental loops and dipoles exhibit figure-eight responses as depicted in Figure 3 
for the azimu th plane. Response maxima are broad but the minima (nulls) are sharp. The 
amplitude response derivative is greatest around the null response position. Therefore, a 
dynamic measure of null position provides the best bearing estimates. 
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Sequential amplitude measurements are made as the directive pattern is rotated in 
a7irmith and the measured bearing is indexed to the azimuth at which some predetermined 
amplitude-null property occurs. 



Figure 3. Azimuth Plane Response Patterns 

The horizontal E-field dipole is rarely used for DF purposes, primarily because most 
signals of interest are vertically polarized and of orthogonal polarization to the horizontal 
dipole [Ref. 6]. 

For comparative amplitude DF, multiple antennas are used to obtain orthogonal 
patterns. Bearing information is obtained from the ratio of the amplitudes of the orthogonal 
patterns. Comparative amplitude DF techniques provide instantaneous bearing information. 

3. Phase Difference 

Bearing measurements using a phase difference require at least two separated 
antennas. A plane wave, arriving at an angle other than the normal to the baseline between 
the two elements, arrives at one element earlier than the other. The time lag between the 
antennas produces a RF phase delay or a differential RF phase between the antenna outputs. 

Figure 4 illustrates the basic phase-delay technique. An incident plane wave arrives 

at an incident angle 6 at antenna 1 inducing a voltage V x exp (jcot). After propagating the 
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distance d sin 9, the incident plane wave induces a voltage V 2 exp (ja> t - & ) in antenna 2, 
where the phase delay given by 

0=(27iJ/X)sin0 . (1) 



Figure 4. Phase-difference DF Parameters 

Therefore the bearing angle is encoded as a function of phase delay. For the phase 
difference technique, the bearing angle is computed by using Eq. 1, where phase difference 
is measured and baseline distance d and wavelength X are known. This phase-delay 
technique can experience phase and bearing ambiguities. If d>/ 1/2 phase ambiguities 
exist; however, bearing ambiguities are resolved by using auxiliary antennas to establish 
known phase references [Ref. 6], 

A variant of the phase measurement method is Doppler direction finding. The 
received frequency from a moving antenna experiences a Doppler shift. If the antenna 
moves along a radial to the emitter position, the Doppler shift is maximized. No Doppler 
shift occurs if the antenna movement is tangential to the direction of propagation. These 
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maxima and minima Doppler shifts may be used to measure the bearing of the incident 


signal. 

4. Time Delay 

An electromagnetic wave propagating between two separated antennas experiences 
a time delay. This time differential of arrival (TDOA) contains bearing angle information. 
In Figure 5, the differential time delay, At, between antennas 1 and 2 is given by 



( 2 ) 


where: 

d = the distance between antennas 1 and 2 (meters); 
c = velocity of light; 

0 = bearing angle (degrees). 

Historically, TDOA has been used for pulsed emissions such as radar. Also, relatively long, 
multiple-wavelength baselines are commonly used [Ref. 6]. 
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IH. ARRAY BEAMFORMING THEORY 


A. OVERVIEW 

Direction finding systems consists of an antenna (or collection of antennas), 
beamforming networks, detectors, and signal processing equipment. It operates by receiving 
and then analyzing the nature of the emitter signal. In noise-limited situations, the source 
signals may consist of a self-product and additive receiver noise. In more complicated 
situations, it may also contain a noise component due to backscatter from the surrounding 
enviro nm ent (clutter). A further source of complication may be the presence of interferences 
produced by jammers beamed at the DF site from directions unknown to operators at the 
site. In the following analysis, unless otherwise noted, these complications are neglected. 
Emitter signals are modeled as time-harmonic plane waves, narrowband assumptions are 
made, and noise is assumed to be white Gaussian. Before examining residue arrays, the 
capabilities of several conventional DF methods are presented. 

B. BASIC SIGNAL AND NOISE MODELS 

Figure 6 shows a uniformly spaced linear array consisting of N elements, and with 
M plane waves arriving at the array from distinct directions. We assume that all these plane 
waves are narrowband with the same carrier frequency/ . The noiseless signal produced at 
the nth element of the array by the wave from source m, denoted as the /nth plane wave, is 
as follows: 


s{n,m,t) = A m cos 


2n f c t + 


N +1 


l”-—J*"* z 


where: 


( 3 ) 
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n = 1,2,...,N (element index) 
m = 1,2... ,M (emitter index) 
f c - carrier frequency 

2 n ? 

k m = vector wavenumber of the mth incident plane wave = — k m 

A 

z = unit vector along the line of the array 
t = time 

The parameter^ denotes the amplitude of the signal s(n,m,t). 



1 2 3 • * * N-l N 

Figure 6. Geometry of an ^Element Linear Array with M Incident Plane Waves 

Let 0 W denote the angle of arrival of the /wth plane wave, measured with respect to 
the no rmal to the array, as illustrated in Figure 6. Then we may express the dot product of 

the two vectors k m and z as follows [Ref. 7]: 
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( 4 ) 


A 2nd . . 

^*^ = — Sm0 m 

where d is the interelement spacing and X is the received radar wavelength. Accordingly, we 
may rewrite Eq. (3) as: 

N + i\ 


s(n,m,t) = 4k cos 


2icfj + 2nj\n- 


2 ) 


sin0. 


( 5 ) 


We may simplify the representation of the signal s(n,m,t) by adopting the complex 
notation for narrowband signals. In particular, we may represent the signal s(n,m,t) by its 
complex amplitude, defined by: 


s(n,m)= A m exp 


f 


u 

2nj l 


N + 1 ' 


n — 


V 


sin0. 


( 6 ) 


Now define the electrical phase angle from element to element along the array as: 


d . 

(j ) m = 27t- sm 0 m . 


( 7 ) 


Also, let a m denote the complex amplitude of the signal s(n,m,t) measured at the center of 
the array; that is: 


(N + \ 




-,m 


( 8 ) 


n=- 


Note that if the number of elements N is odd, the center of the array represented by 
N + — corresponds to the center element. On the other hand, when N is even, the center 


of the array is a fictitious point. In any event, we may use Eq. 7 to rewrite Eq. 6 in the 
simplified form 


s(n,rn) = a m exp 


j K 


N+i 
2 J 




( 9 ) 


The signal measured at the output of any element in the array differs from the signal 
actually received by that element by an amount attributed to noise. In practice, this noise is 
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modeled as a white ergodic random process. To limit the effects of the noise, it is customary 
to filter the output of each element in the array to the same narrow band of frequencies 
occupied by the actual received signal. As a result of this operation, we have an observed 
cignal which for each element of the array, consists of the actual received signal plus a 
narrowband noise. For element n we may describe this narrowband noise by 

w(«,0 = B„cos(2n/ c f + |5„) (10) 

where the amplitude B n is a Rayleigh distributed random variable and the phase is 
uniforml y distributed random variable over the range (0, 2%). Representing the narrowband 

noise w(n,t) by its complex amplitude, w{n ), we have: 

w(n) = 5„exp(yp„) . (11) 

This is a complex-valued random variable that is Gaussian-distributed with a mean that is 


typically zero. 

Equation 9 defines the complex amplitude of the signal received by element n due 
to the mth. plane wave. Summing the contributions produced by a set of M plane waves, and 
then adding to the result the complex amplitude of the narrowband noise, we may express 
the complex amplitude of the observed signal as: 

;c(n) = s(«) + w(n) 

M 

= Y 4 s(n,m) + w(n) (12 ) 

m= 1 
M 

=2X ex p 

m= 1 

We may put Eq. 12 into a more compact form by using matrix notation. First introduce the 
following definitions: 


J\ *■ 


N+i 
2 J 


♦. 


+w(«), n = l,2,...,N 


1. AnNxJ vector 
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x= 


called the observed signal vector. 


2. An N x 1 vector 


*(D 

4A 


x(l) 

x(2) 

x(N) 


S = 


4n) 


2>M 

m= 1 
M 

Z 5 ( 2 ’ m ) 


m -1 
M 


Y,s( N > m ) 


,m= 1 


called the received signal vector. 


3. An Mxl vector 


A = 


a(l) ' 

a(2) 

a(M) 


called the signal-in-space vector. 


4. An jVxM matrix 


(13) 


(14) 


(15) 


exp (jkd cos 0 ,) 


exp (jkd cos 0 2 ) 


exp (jkd cos 0 M ) 


exp (jk(N-l)dcos0 ,) exp (jk(N-l)dcos0 2 ) ••• exp (y'£(A-l)<icos<9 M ) 

(16) 

of propagation delays called the direction matrix. 
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5. An N x 1 vector 


f w (0 



w(N) 


called the noise vector. 

Then, using these definitions, we may rewrite Eq. 12 in the form 

x = s+w 

or, equivalently, 

X = BA + W 


(17) 


(18) 

(19) 


where 

S = BA . (20) 

The observed signal vector X represents a snapshot of data corresponding to a particular 
instant of time. Ordinarily, several independent measurements are made, so that the data 
available for processing may be expressed in the form: 

X(i) = BA(/) + W(i), i = 1,2,...,/ (21) 

where / is the total number of snapshots taken. 

Based on such a set of data, we may use temporal averaging to improve the 
estimation of the angles of arrival of the plane waves incident upon the array. The signal-in¬ 
space vector A(/) is modeled as a stochastic process because the behavior of the sources 
responsible for its generation is, in general, unpredictable. 

The spatial correlation matrix of the observed signal vector X is defined by 

R = E[x , (i)X r (>)] m 

where E is the expectation operator, the asterisk denotes complex conjugation, and the 
superscript T denotes transposition. Substituting Eq. 21 in Eq. 22, and recognizing that the 
signal-in-space vector A and the noise vector W are statistically independent, we get 
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( 23 ) 


R = E[B*A*A r B r ] + E[w*W r ] 

= B*R a B r +R W 

where: 

R a = spatial correlation matrix of the signal - in - space vector A 
= E[A*A r ] 

R w = spatial correlation matrix of the noise vector W 
= E[W*W r ] . 

The matrix product 

R s = B*R a B r (24) 

represents the spatial correlation matrix of the signal component of the signal vector S. 

With the noise vector W assumed white noise, we have 

R W = ( 25 ) 

where o\ is the variance of the elemental noise w(n) for all n, and I represents the identity 
matrix. Depending on the structure of the spatial correlation matrix R a , we may distinguish 

two different cases: 

1. Jointly Uncorrelated 

The sources responsible for the signal-in-space vector A are jointly uncorrelated or 
noncoherent so that the off-diagonal elements of the corresponding spatial correlation 
mat rix R a are all zero. Then we find that the spatial correlation matrix of the received signal 
vector, R s = B* RaB f is a Toeplitz matrix. A square matrix is Toeplitz if all the elements on 
its main diagonal are equal, and likewise for the elements on any other diagonal parallel to 
the main diagonal. The received signal vector S is spatially stationary if its spatial 
correlation matrix is Toeplitz, or vice versa. Clearly, if the spatial correlation matrix R s is 
Toeplitz, the spatial correlation matrix of the observed signal vector, R is likewise Toeplitz. 
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2. Correlated Coherent 

The sources responsible for the received signal vector S are correlated coherent. The 
result is that the spatial correlation matrix R s is non-Toeplitz, and the received signal vector 
S is said to be spatially nonstationary. Clearly, If R is non-Toeplitz, the spatial correlation 
matrix, R, pertaining to the observed signal vector, is also non-Toeplitz. 

C. ARRAY MATRIX DECOMPOSITION FOR MULTIPLE EMITTER DOA 


For the location of the directions of arrival, the autocorrelation matrix takes the form 
of a Toeplitz matrix. 

X = BA + W 

A = [a(l),tf(2),...,a(M)f (26) 

W = [ w(l), w{2),w( N)f 


Vector X is explained in Eq. 13 and vector B is defined in Eq. 16; it contains the 
propagation delays, that is B = [U / ,U 2 ,...,U A/ ]where U m =V{0 m ), and, for equally 
spaced elements. 



~ e jM cose 


l 

gjkd cos 0 

jkd N cos 6 


gjk(N-\)dcos9 


(27) 


Thus, the signal correlation matrix can be written as 

r = e[b*a*a'b']+e[w*w'] 
= B*R a B'+R w 


(23) 


where t stands for the complex conjugate transpose. For an equally spaced array with N 
elements located at 0, d, 2d,...,{N-l)d the size of the correlation matrix in Eq. 23 is NxN. In 
this case, the maxim um number of emitters that can be detected cannot exceed (N-l). 

The matrix R is Hermitian and always positive definite if thermal element noise is 
present Thus R can be diagonalized using a rotation matrix E where 
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ERE' = 


(28) 


A,, ••• 0 

• • • 

• • • 

0 ••• X N _ 

The matrix E is unitary (EE'=I) and its columns are the eigenvectors of R, where 

The signal covariance matrix is diagonalized by the same rotation matrix. For example, if 
there is only one signal present (N=l), the weight vector for the single signal case is: 

G = (constant)e,. 

A measure of the array response is given by the spectral estimator 


m= 


i 

|G'u(e )| 2 



(30) 


A large value of P denotes high array sensitivity; that is a pattern notch or null. Note 
that this is not the same condition that provides good signal to noise for a conventional 
array. If more than one emitter is present the spectral estimator becomes 


1 

- 7 - a ( 31 ) 

[ e n e 2»"’» e A'-Af] u(0)| 

The power response for the conventional equally spaced array of 6 elements using ordinary 
beamforming is shown in Figure 7. In this case there is no ability to resolve two signals, the 
null locations are fixed relative to each other. Figure 8 shows the response of a 
Caratheodory array [Ref. 1] with four elements located at 0, d, 5d/2, and 3d. The emitter 
directions are at 20, 50, 90, 130 and 160 degrees respectively. Note that all of the emitters 
are essentially resolved, but the sensitivity is only lOdB in some cases. 
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Figure 7. Conventional Equally Spaced Array with Six Elements when the Emitter 

Directions are 20,50,90,130, and 160 
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Figure 8. Conventional Processing Scheme for a Caratheodory-array with Four 
Elements when the Emitter Directions are 20,50,90,130, and 160 
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IV. RESIDUE ANTENNA ARCHITECTURE 


A. OVERVIEW OF RESIDUE NUMBER SYSTEMS 

In the first century A.D. the Chinese scholar Sun Tzu authored a book posing the 
problem of determining a number having remainder 2, 3, and 3 when divided by 3,4, and 5 
respectively. The answer is 23 and is found by using a three-moduli (3,4,5) number system. 
Later this theorem was called the Chinese Remainder Theorem [Ref. 8]. Eventually this 
theorem became an important cornerstone in the modem theory of residue number system 
(RNS) arithmetic. 

Digital systems that use residue arithmetic units play an important role in ultra-high 
speed, dedicated, real-time systems that support pure parallel processing of integer valued 
da ta Residue arithmetic performs addition, subtraction, and multiplication as concurrent 
operations, side-stepping one of the principal arithmetic delays: managing carry 
information. Some applications such as the fast Fourier transform, (iigital filtering, and 
image processing utilize the efficiencies of RNS arithmetic in addition and multiplication. 

The RNS has been employed efficiently in the implementation of digital processors. 
The RNS is composed of a number of moduli m v The integers within each RNS modulus 
are representative of a sawtooth waveform with the period of the waveform equal to the 
modulus. 

Figure 9 shows the RNS folding waveform and the corresponding RNS states for 
the moduli /w,=3, m 2 = 4, m 3 = 5. In this scheme, given m u the integer values within the 
individual modulus are given by the row vector: 

[0,l,...,w ( -l] (32) 

The required number of comparators for the zth channel is m— 1. Notice that all the integers 
within Eq. 32 form a complete system of length m l . The dynamic range of this scheme, as 
well as the number of quantization levels, is given by: 
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( 33 ) 


N 

Dynamic Range = n 

(=i 


Table 1 illustrates the RNS encoding values for the moduli (3,4,5). 



m l =3 


m 2 =4 


m 3 = 5 


Figure 9. RNS Folding Waveforms and Output Code for the Moduli (3,4,5) 
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Analog Moduli 

Input *w,=3 0*2=4 *n 3 =5 

0 

0 

0 

0 

1 

1 

1 

1 

2 

2 

2 

2 

3 

0 

3 

3 

4 

1 

0 

4 

5 

2 

1 

0 

6 

0 

2 

1 

7 

1 

3 

2 

8 

2 

0 

3 

9 

0 

1 

4 

10 

i I 

2 

0 

11 

2 

3 

1 

12 

0 

0 

2 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

47 

2 

3 

2 

48 

0 

0 

3 

49 

1 

1 

4 

50 

2 

2 

0 

51 

0 

3 

1 

52 

1 

0 

2 

53 

2 

1 

3 

54 

0 

2 

4 

55 

1 

3 

0 

56 

2 

0 

1 

57 

0 

1 

2 

58 

1 

2 

3 

59 

2 

3 

4 

60 

0 

0 

0 


Table 1. RNS Encoding Procedure for the Moduli (3,4,5) 


B. RNS ANTENNA ARCHITECTURE 

In the RNS antenna architecture shown in Figure 10, phase detectors are used to 
combine the individual pairs of elements. The comparator threshold levels and subsequent 
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logic block are particular to the RNS. The phase response (A^) of each interferometer 
corresponds to a sawtooth folding waveform. 



Figure 10. RNS Antenna Architecture 
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For a given modulus m p the number of folds needed to apply the RNS code within 
the dynamic range is given by: 

dynamic range M 


W i,RNS - 


m : 


m ; 


(34) 


From Figure 10 


9 = 0°=> A<f> = 0° , 

9 = 90°=> A^ = — 


(35) 


9 = -90° => = - 


2nd 

~Y 


For a DOA of 180°, the phase response corresponds to a change in phase of 2n. We require 
that 

A^90°) - A^-90°) = 2 n , (36) 

which implies 


4tt d _ , ^ 

-= 2tc, or a = — . 

X 2 

Finally, the distance between elements is given by: 

Ck\ M(V 

“/,RNS - n i, RNS 


(37) 


(38) 


\2) 

Using the same moduli set with center frequency of f 0 = 8.5 GHz, the distances for each 


pair of elements are: 


4=m , 
d 2 = 7.5X , 
d 2 = 6.0A, . 


These interferometer distances are used to investigate numerically the performance of the 
RNS antenna architecture. 


27 


C. RNS ARRAY CORRELATION MATRIX FOR MULTIPLE EMITTERS 

In this section we briefly investigate the capability of the RNS array with m,=3, 
m 2 =A, and m 3 =5 to resolve multiple emitters at the same frequency. The same approach is 
used for as the minimum redundancy array (Chapter 2, Section C). For a minimum 
redundancy array, not all elements must be present, but all correlation lags are present (d, 
2d, 3d„ .... (N-l)d) since all spacings are present at least once. Therefore the array correlation 

matrix can be filled and consequently up to —— emitters can be resolved with N 

elements. 

For the RNS array not all elements spacings are present. In the case of m x =3, m 2 =4 
and mf=5 the spacings that are present are 0 d, 12 d (6A), 15 d (7.5 A), 20 d (10A), 3 d 
(1.5 A), 5 d (2.5 A ), and %d (4 A) as shown in Figure 11. 


l(U(20d) 


-7.5A(15d) 

62(12d)- 


^ V V V 


r 1 

s \ 

S ' 


<1.5A> 

<- 

<—2.5/1-> 

-4A -> 


Figure 11. RNS Antenna Spacing for the (3,4,5) Array 
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The standard periodic array would have a completely filled 21 x 21 correlation 
matrix. Not all spacings are represented in the RNS array, and if the missing correlation lags 
are set to zero (or noise) the matrix would be: 


R = 


Ro 0 0 R 3 0 R 5 0 0 R g 0 

0 Ro 0 0 R 3 0 R 5 0 0 R 8 

0 0 R 0 0 0 R 3 0 R 5 0 0 

R’ 3 0 0 Ro 0 0 R 3 0 R 5 0 

0 R* 3 0 0 Ro 0 0 R 3 0 R 5 

R’s 0 R’ 3 0 0 Ro 0 0 R 3 0 

0 R% 0 R% 0 0 Ro 0 0 R 3 

0 0 R* 5 0 R*3 0 0 Ro 0 0 

R* g 0 0 R* 5 0 R* 3 0 0 Ro 0 

0 R’g 0 0 R's 0 RY 0 0 Ro 

0 0 R’g 0 0 R* 5 0 R*3 0 0 

0 0 0 R’g 0 0 R* 5 0 R* 3 0 

r* 12 o o o r* 8 o o r* 5 o r* 3 

0 RY 0 0 0 R’g 0 0 R* 5 0 

0 0 RY 0 0 0 R’g 0 0 R*5 

RY 0 0 RY 0 0 0 R’g 0 0 

0 RY o 0 RY o 0 0 R’g 0 

0 0 RY o 0 RY 0 0 0 R’g 

0 0 0 RY o 0 RY 0 0 0 

0 0 0 0 RY 0 0 RY o 0 

RY 0 0 0 0 RY 0 0 RY 0 


0 0 R 12 0 0 R,s 0 0 0 

0 0 0 R 12 0 0 R,5 0 0 

Rg 0 0 0 R, 2 0 0 R 15 0 

0 Rg 0 0 0 R }2 ^ o R-i 5 

0 0 Rg 0 0 0 R 12 0 0 

R 5 0 0 Rg 0 0 0 R 12 0 

0 R 5 0 0 Rg 0 0 0 R 12 

R 3 0 R 5 0 0 Rg 0 0 0 

0 R 3 0 R 5 0 0 Rg 0 0 

0 0 R 3 0 R 5 0 0 Rg 0 

Ro 0 0 R 3 0 R 5 0 0 Rg 

ORoO 0R 3 0R 5 0 0 

0 ORoO 0R 3 0R 5 0 
R’ 3 0 0 Ro 0 0 R 3 0 R 5 

0 R’ 3 0 0 Ro 0 0 R 3 0 

R's 0 R’ 3 0 0 Ro 0 0 R 3 

0 R’s 0 R’ 3 0 0 Ro 0 0 

0 0 R’ 5 0 R’ 3 0 0 Ro 0 

R’s 0 0 R’s 0 R’ 3 0 0 Ro 

0 R’g 0 0 R% 0 R* 3 0 0 

0 0 R’g 0 0 R’s 0 R’ 3 0 


0 R 2 o 
0 0 
0 0 
0 0 
R,5 0 
0 R,s 
0 0 
r 12 0 
o r 12 
0 0 
0 0 
r 8 o 

0 Rg 
0 0 
Rs 0 
0 R 5 
R 3 0 
0 R 3 
0 0 
Ro 0 
0 R 0 


Figure 12 shows the array response for 50, 90 and 130 degrees. As expected the 
performance is not good because of the large number of zeros in the array. 
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Figure 12. RNS Array Response for Multiple Emitters when the Emitter Directions 

are 50,90, and 130 

D. COMPUTER SIMULATION OF RNS CODE 

The quantization of the phase response from each pair of elements is made separately using 
a small bank of comparators. The threshold levels for each comparator are shown in Table 
2. If the phase response is plotted as a function of sin0, an ideal sawtooth waveform is 
obtained as shown at the top of Figure 13. An actual plot of the phase response is shown in 
Figure 14 as a function of the input direction of arrival 0. The sawtooth waveform stretches 
out as 0, the input direction of arrival, goes from 0 to 90°. For a specific input direction 
of arrival, there will be a corresponding integer number of comparators in the “on” 
state within the modulus tw,-. The value of the integers depends on the measured phase 
difference associated with each pair of elements for the specific input direction of arrival. 
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The three RNS codes (one for each pair of antenna elements) are indices to an RNS logic 
look-up table that gives the resolved direction of arrival. Table 3 shows the three RNS 
codes and the resolved direction of arrival. The numbered codes represent the number of 
comparators in the “on” state. 


Moduli 

T, 

t 2 

T, 

t 4 

3 

-1.047 

1.047 



4 

-1.571 

0.00 

1.571 


5 

-1.885 

-0.628 

0.628 

1.885 


Table 2. Threshold Levels for the Moduli (3,4,5) 



Input Direction of Arrival (Sin 0) 

Figure 13. Phase Response as a Function of sin0 for the Moduli (3,4,5) 
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Figure 14. Phase Response as a Function of the Input Direction of Arrival (0) 

for the Moduli (3,4,5) 

Figure 15, shows a plot of the resolved direction of arrival as a function of 0, the 
input direction of arrival. The RNS antenna system gives the direction of arrival 
information over the complete field of view without ambiguity (DOA increment is 0.3°). 
However the accuracy varies depending on the actual direction of arrival of the incoming 
signal. Notice also that the quantization width increases as the direction of arrival 
approaches ±90°. 
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Resolved Direction of Arrival (degrees) 



Figure 15. Resolved Direction of Arrival in Function of the Input Direr 

when DOA Increment is 0.3° 
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V. ENCODING ERROR CORRECTION 


A. ENCODING ERRORS 

The DO A estimates previously shown in Figure 15 are based on ideal hardware 
performance. One significant assumption is that all of the comparator thresholds are ideal 
and cross the falling phase response waveforms together at each DOA code transition point. 
If the samples are generated by stepping the emitter direction between -90° and 90° in 
increments of 0.03°, for example, some samples will fall about the code transition point and 
result in encoding errors due to some comparators do not switch simultaneously. When the 
emitter is swept through an angular sector, the DOA estimates from a realizable system 
could contain large spikes or “glitches” at these points. An illustration of the encoding 
errors is given in Table 4. The columns in Table 4 correspond to the first 5 RNS codes 
shown in Table 3 with each code sampled three times. The seventh sample shows the nature 
of an encoding error where some of the comparators in modulus 5 do not switch at the same 
time. The seventh sample should be m x = 2, m 2 =2 and m 3 =2. This results in a glitch at this 
sample within the transfer function as shown in Figure 16. 


Sample 

1 

m 

rn 

rn 

5 

6 


B 


10 

11 

EEI 

13 

14 

Modulus 3 

0 

o 

o 

B 

1 

1 

2 

B 

B 

ra 

0 

D 

0 

1 

Modulus 4 

D 

0 

o 

D 

1 

1 

2 

B 

2 

B 

3 

B 

0 

0 

Modulus 5 

JL 

O 

o 

_L 

B 

1 

1 

B 

2 

B 

3 


■ 4 

4 


Table 4. Encoding Error in the RNS Output Code 


To examine the nature of the encoding errors, a simulation was performed with the 
threshold levels shifted from the ideal (error-free) threshold levels in Table 2 to a set of new 
threshold levels. The errors are induced in order to compare three interpolation methods 
designed to eliminate the encoding errors. The new threshold values are shown in Table 5, 
and the glitches in the DOA are shown in Figure 17, with the total number of glitches equal 
to 44. 
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Figure 16. RNS-to-Decimal Logic Block and Glitch 
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Moduli 

T, 

t 2 

t 3 

t 4 

3 

-1.047 

1.047 



4 

-1.560 

0.000 

1.560 


5 

-1.870 

-0.640 

0.640 

1.870 


Table 5. New threshold Levels for the Moduli (3,4,5) 


Encoding Error in Quatization (glitch= 44,sampling distanca:ds=0.03 degree) 
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Figure 17. Glitches and Resolved Direction of Arrival 


B. ADDITIONAL COMPARATORS TO ISOLATE ENCODING ERRORS 

The original threshold levels within the minimum modulus are set to the level of 
1.047, -1.047 radians. The phase response for the minimum modulus (m,= 3) is shown in 
Figure 18 and illustrates the position of the comparator threshold levels. 
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Figure 18. Original Threshold Level in Modulus 3 

To isolate the possible encoding errors, 2m, comparators are used in the minimum 
modulus (instead of m, -1). The comparator thresholds levels are such that a small band is 
set up about the code transition points. The values of the threshold levels depend on the 
width of the band and can be expressed as a percentage of the least significant angle (LSA) 
width. The six-threshold levels for modulus m,=3 are shown in Figure 19, and listed in 
Table 6 as a function of the LSA percentage. 

The encoding errors are isolated by comparing the parity value for the 2m, 
comparators in the minimum modulus. If the parity value is odd number (1, 3 and 5), then 
the sample falls outside the small transition band and does not present a problem. If the 
parity is an even number (2,4 and 6), the sample represents a possible DOA encoding error. 
The goal of the interpolation method is to determine the best DOA estimate given the fact 
that the sample could present a possible error. 
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Plase Output 



Figure 19. Additional Threshold Levels and RNS Output Code 


LSA Percentage 

T„ 

T,* 

T* 

T a4 

T,s 

T, 6 

10.00 

-3.0369 

-1.1519 

-0.9425 

0.9425 

1.1519 

3.0369 

9.00 

-3.0473 

-1.1414 

-0.9529 

0.9529 

1.1414 

3.0473 

8.00 

-3.0578 

-1.1310 

-0.9634 

0.9634 

1.1310 

3.0578 

7.00 

-3.0683 

-1.1205 

-0.9739 

0.9739 

1.1205 

3.0683 

6.00 

-3.0788 

-1.1100 

-0.9844 

0.9844 

1.1100 

3.0788 

5.00 

-3.0892 

-1.0996 

-0.9948 

0.9948 

1.0996 

3.0892 

4.00 

-3.0997 

-1.0891 

-1.0053 

1.0053 

1.0891 

3.0997 

3.00 

-3.1102 

-1.0786 



1.0786 

3.1102 

2.00 

-3.1206 

-1.0681 

-1.0263 

1.0263 

1.0681 

3.1206 

1.00 

-3.1311 

-1.0577 

-1.0367 

1.0367 

1.0577 

3.1311 


Table 6. LSA Percentage and Six-Additional Threshold Values 
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C. INTERPOLATION 


A block diagram of the RNS antenna architecture including the interpolation 
processing is shown in Figure 20. 
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Figure 20.The Procedure of DOA Decoder 

The RNS-to-binary logic block is set after the parity circuit as the interpolation 
procedures are accomplished in the binary number system as described below. The RNS 
code and corresponding binary values are shown in Table 7. The encoding errors at samples 
7 and 13, and the RNS code in Table 4 are converted to a six-bit binary code and are shown 
in Table 8. 
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Angle 

RNS 

Code 

Six-Bit 

Binary 

Code 

Angle 

RNS 

Code 

Six-Bit 

Binary 

Code 

-82.5824 

0 

0 

0 

000000 

0.9551 

0 

2 

0 

011110 

-72.0627 

1 

1 

1 

000001 

2.8664 

1 

3 

1 

011111 

-66.5593 

2 

2 

2 

000010 

4.7809 

2 

0 

2 

100000 

-62.1158 

0 

3 

3 

000011 

6.7007 

0 

1 

3 

100001 

-58.2581 

1 

0 

4 

000100 

8.6382 

1 

2 

4 

100010 

-54.7864 

2 

1 

0 

000101 

10.5655 

2 

3 

0 

100011 

-51.5928 

0 

2 

1 

000110 

12.5152 

0 

0 

1 

100100 

-48.6111 

1 

3 

2 

000111 

14.4797 

1 

1 

2 

100101 

-45.7968 

2 

0 

3 

001000 

16.4618 

2 

2 

3 

100110 

-43.1187 

0 

1 

4 

001001 

18.4644 

0 

3 

4 

100111 

-40.5534 

1 

2 

0 

001010 

20.4907 

l 

0 

0 

101000 

-38.0832 

2 

3 

1 

001011 

22.5442 

2 

1 

1 

101001 

-35.6940 

0 

0 

2 

001100 

24.6287 

0 

2 

2 

101010 

-33.3745 

1 

1 

3 

001101 

26.7487 

i 

3 

3 

101011 

-31.1155 

2 

2 

4 

001110 

28.9091 

2 

0 

4 

101100 

-28.9091 

0 

3 

0 

001111 

31.1155 

0 

1 

0 

101101 

-26.7487 

1 

0 

1 

010000 

33.3745 

1 

2 

1 

101110 

-24.6287 

2 

1 

2 

010001 

35.6940 

2 

3 

2 

101111 

-22.5442 

0 

2 

3 

010010 

38.0832 

0 

0 

3 

110000 

-20.4907 

1 

3 

4 

010011 

40.5534 

1 

1 

4 

110001 

-18.4644 

2 

0 

0 

010100 

43.1187 

2 

2 

0 

110010 

-16.4618 

0 

1 

1 

010101 

45.7968 

0 

3 

1 

110011 

-14.4797 

1 

2 

2 

010110 

48.6111 

1 

0 

2 

110100 

-12.5152 

2 

3 

3 

oioili 

51.5928 

2 

1 

3 

110101 

-10.5655 

0 

0 

4 

011000 

54.7864 

0 

2 

4 

110110 

-8.6382 

1 

1 

0 

011001 

58.2581 

1 

3 

0 

110111 

-6.7007 

2 

2 

1 

BTITiTTilM^mgai 

IB 

0 

1 

111000 

-4.7809 

0 

3 

2 

011011 

66.5593 

0 

1 

2 

111001 

-2.8664 

1 

0 

3 

011100 

72.0627 

1 

2 

3 

111010 

-0.9551 

2 

1 

4 

011101 

82.5824 

IEKKI 

111011 


Table 7. RNS code and Six-bit Binary Code in Moduli (3,4,5) 
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Sample 

1 

2 

3 

4 

5 

6 

|7 

8 

9 

10 

[]L 

E9 

13 

14 

MSB 
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0 

0 

0 

0 

0 

o 

0 

0 

ra 

O 

0 

0 

- 

0 

o 

0 

0 

D 
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ra 

o 

O 

o 

O 

1 

0 

- 

0 

0 

0 

D 

0 

0 

0 

El 

0 

O 

El 

o 

1 

0 

- 

0 

0 

O 

o 

0 

0 

1 

O 

0 

0 

El 

ra 

0 

1 

* 

o 
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0 

O 

O 

ra 


UL 

1 

1 

1 

n 

0 

0 

LSB 

D 

D 

O 

D2 

m 

□3 

mm 

El 

0 

1 

1 

LU 

0 

0 


MSB represents the Most Significant Bit of the converter 


LSB is the Least Significant Bit in the six-bit converter 

Table 8. Encoding Error in Six-bit Binary Code 

1. LSB Shift Method 

In Table 4 glitches occured at samples 7 and 13. From that binary data the output 
signal jumps from the value given by 000001 through 000110 to 000010 as shown in Table 
8. To resolve these glitches, we investigate the use of a 7x5 Error Correction Block (ECB) 
matrix. The output of a parity circuit determines whether the sample is good or bad. If the 
parity output is an odd number (1, 3 and 5) the sample is good, whereas an even number (2, 
4 and 6) the output represents a possible bad sample. The parity circuit and error correction 
block are shown in Figure 21. 

a. LSB Shift Algorithm 

The parity and the 6-bit word (from the PLA) are clocked in parallel through 
a 5-state buffer. The buffer allows the 6-bit word to be replaced with an interpolated value 
depending on the value of the parity. If the parity is odd (good sample), the 6-bit word is 
clocked through the buffer unaltered. If the parity is even, the values ahead and behind the 
even parity value are used to interpolate the possible encoding error. The interpolation 
process compares the most significant bits (MSBs) of the odd parity words stopping at the 
point where the MSBs do not agree. The MSBs that are the same are copied into the MSBs 
for the even parity word. The LSBs come directly from the original (even parity) word. 
Figure 22 shows an example where there is one out of five possible bad samples. In the 
interpolation process, P 2 and P 0 combine to form the MSBs of T, (an intermediate state). 
The LSBs of T, come from Pj. At time t+1, T, is shifted into P 2 . The situation where two 
samples are even parity is shown in Figure 23. In this case P 3 and P 0 combine to form the 
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MSBs of Tj and T 2 . The LSBs of T, and T 2 are copied directly from P, and P 2 . T, is shifted 
into P 2 and T 2 is shifted into P 3 at t+1. The interpolation procedure for three out of five 
samples having even parity is shown in Figure 24. The interpolation procedure for four out 
of five samples having even parity is shown in Figure 25. 


Parity = 1 (Sample OK) 
Output 0 (Sample Bad) 


7x5 Error Correction Block 


> 


LSB _ 

Sn= Shift Register 



Figure 21. Parity Circuit and Error Correction Block 
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b. Simulation Results 

The DOA transfer function including the encoding error? was shown in 
Figure 17. The LSB shift method with the interpolation width equal to 1% of the LSA is 
shown in Figure 26. The total number of glitches has been reduced from 44 to 42. Figure 27 
through Figure 30 show the result of the LSB shift method with the interpolation width 
ranging from 2% to 5% respectively. Figure 31 summarizes the results and shows the total 
number of glitches as a function of the isolation bandwidth (percentage of LSA). The larger 
the bandwidth, the fewer the number of glitches. 


Interpolation Result, 1 % LSB-Shift Method (glitch=42, sampling distance:ds=0.03 degree) 



-80 -60 -40 -20 0 20 40 60 80 

Input Direction of Arrival (degrees) 


Figure 26. LSB-Shift Interpolation Results, 1% of LSA 
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Figure 28. LSB-Shift Interpolation Results, 3% of LSA 
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Resolved DOA and Quantization Error (degrees) Resolved DOA and Quantization Error (degrees ) 


Interpolation Result,4 % LSB-Shift Method (glitch= 15, sampling distance:ds=0.03 degree 



Figure 29. LSB-Shift Interpolation Results, 4% of LSA 


Interpolation Result.5 % LSB-Shift Method (glitch= 8, sampling distance:ds=0.03 degree) 



Figure 30. LSB-Shift Interpolation Results, 5% of LSA 
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Figure 31. Number of Glitches vs LSA Percentage After LSB-Shift Method 

2. Random LSB Method 

The random LSB method compares the most significant bits (MSBs) of the odd 
parity words stopping at the point where the MSBs do not agree. The MSBs that are the 
same are copied into the MSBs for the even parity word. The LSBs do not come from the 
original (even parity) word. Instead the LSBs are chosen randomly (0 or 1). The 
interpolation procedure for three out of five samples having even parity is shown in Figure 
32. 

The simulation for the random LSB method with the interpolation width changing 
from 1% to 5% are shown in Figure 33 to Figure 37. The total number of glitches changed 
randomly, so it is difficult to determine any relationship between the LSA percentage and 
glitches. The results are summarized in Figure 38 and the randomizing of the LSBs seems 
to create more problems. 
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Interpolation Result.l % Random-LSB Method (glitch- 82, sampling distance:ds=0.03 degrees 



Interpolation Results % Random-LSB Method (glitch= 95, sampling distance:ds=0.03 degrees 



Figure 34. Random-LSB Method Result, 2% of LSA 
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Resolved DOA and Quantization Error (degrees) 


Interpolation Result,3 % Random-LSB Method (glitch= 114, sampling distance:ds=0.03 degree) 



Figure 35. Random-LSB Method Result, 3% of LSA 

Interpolation Result,4 % Random-LSB Method! glitch* 98. sampling distance:ds=0.03 degree) 



Figure 36. Random-LSB Method Result, 4% of LSA 
























































Number of Gliches 


Interpolation Result.5 % Rendom-LSB Method (giitch= 67, sampling distance:ds=0.03 degree) 



Figure 37. Random-LSB Method Result, 5% of LSA 



Figure 38. Number of Glitches vs LSA Percentage after Random LSB Method 


























3. Shift Last Good-Sample Method 

In the last two methods, some error samples cannot be resolved completely. This is 
because some portion of least significant bits remains in the error correction block buffer 
without changing. With the previous methods some glitches may remain in the sample. 
However, the shift last good-sample method removes the encoding error perfectly after 
interpolation by shifting directly from good sample's 6-bit value to error sample's 6-bit 
value. The procedure of this method is shown in Figure 39. 


56 



Figure 39. Procedure of Shift Last Good-Sample Method 


























Interpolation Resuft.l % Shift Last-Good Bit Method (glitch= 36, sampling distance: ds=0.03 degrees) 



Figure 40. Shift Last Good-Sample Method Result, 1% of LSA 


Interpolation Result.2 % Shift Last-Good Bit Method (glitch= 23, sampling distance: ds=0.03 degrees) 



Figure 41. Shift Last Good-Sample Method Result, 2% of LSA 
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Interpolation Results % Shift Last-Good Bit Method (g!itch= 9, sampling distance:ds=0.03 degrees) 



Figure 42. Shift Last Good-Sample Method Result, 3% of LSA 

Interpolation ResuM % Shift Last-Good Bit Method (glitch= 7, sampling distance: ds= 0.03 degrees) 



Figure 43. Shift Last Good-Sample Method Result, 4% of LSA 
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Interpolation Results % Shift Last-Good Bit Method (glitch= 3, sampling distance:ds=0.03 degrees) 



Figure 44. Shift Last Good-Sample Method Result, 5% of LSA 



Percentage of LSA 


Figure 45. Number of Glitches vs LSA Percentage After Shift Last Good-Sample 

Method 
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VI. CONCLUDING REMARKS 


The important contribution of this thesis is detailing a residue antenna architecture 
and interpolation methods that may hold the key to resolving the encoding errors (glitches) 
using an error correction block. The essential part of interpolation procedure is to determine 
the additional threshold levels according to LSA percentage in the first modulus (1^=3), so 
that the parity circuit can be applied to the error correction logic block system. Using the 
residue number system and the six-bit binary code system it is very fast and develops good 
resolution for converting the analog input into the direction of arrival (DO A) output. 

The three interpolation method results show that the number of glitches is 
determined by the percentage of LSA in three trials. In the LSB shift method and shift last 
good sample method, the simulation results indicate that the higher the LSA percentage the 
better the result (i.e. fewer glitches after interpolation). The random LSB method shows no 
trends between LSA percentage and the interpolation results because of the number of 
glitches is randomly changed. A comparison of three interpolation methods is shown in 
Figure 46. 

The simulated results indicate that the RNS provides a more robust solution to high- 
resolution DF and also has excellent performance. In fact, since the phase response has a 
"uniform" distribution as a function of angle of arrival it is possible to obtain the direction 
of arrival information over the entire field of view. To avoid ambiguities in the resolved 
direction of arrival the spacing between elements must be established for the highest 
frequency emitter. On the other hand, the frequency of the incoming signal has to be known 
in order to make the necessary corrections in the RNS logic block using a fast correction 
algorithm to obtain the correct direction of arrival. The performance of the algorithm was 
quantified and shown to have excellent results over a large bandwidth. Future work should 
include a demonstration array with both the microwave and A/D hardware incorporated as 
shown in Figure 10. 
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Percentage of LSA 

Figure 46. Comparison of the Three Interpolation Method Results 









APPENDIX MATLAB CODE 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
% This program 'Phase3' calculates the Direction Of Arrival (DOA) using % 

% the phase response of the array and RNS code in fast A/D converter % 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


% program 1 
clear all 

% Calculate the Dynamic Range 

ml=3; 

m2=4; 

m3=5; 

M=ml*m2*m3; 

rad=pi/180; 

ds=3\input('What is the step size (degree) ?') 

LSA_percent ==input('What is the percentage of LSA ?(%)') 
fo=input('What was the frequency of the antenna (GHz) ?') 
fr=input('What is the frequency of the received signal(GHz)?') 
fo=fo*le9; 
fi=fr*le9; 

theta=linspace(-90,90,M/ds); 
thet=sin(theta*rad)* 90; 
c=3e8; 

L=c/fo; 

Ll=c/fr; 

% Distances Between Each Pair of Elements 

dl=0.5*M/ml*L; 

d2=0.5*M/m2*L; 

d3=0.5*M/m3*L; 

% Generates the Phase Response of Each Pair of Elements 

phi 1 =(2*pi*dl/L 1 )* sin(theta*rad); 

mod 1 p=angle(exp(j * ((phi 1 +pi)))); 

phi2=(2*pi*d2/Ll)*sin(theta*rad); 

mod2p=angle(exp(j *((phi2)))); 

phi3=(2*pi*d3/L 1 )*sin(theta*rad); 

mod3p=angle(exp(j*((phi3)+pi))); 

% Find the Threshold Values for Each Modulus 
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thresl 


% Generates the Comparators Output RNS code for Modulus m 1 =3 ,m2=4,m3=4 

for i=l:M/ds 

if modlp(i)<=Tl(l) 

modl(i)=0; 

elseif modlp(i)>Tl(2) 

modl(i)=2; 

else 

modl(i)=l; 

end 

if mod2p(i)<=T2(l) 
mod2(i)=0; 

elseif mod2p(i)<=T2(2) 
mod2(i)=l; 

elseif mod2p(i)<=T2(3) 

mod2(i)=2; 

else 

mod2(i)=3; 

end 

if mod3p(i)<=T3(l) 
mod3(i)=0; 

elseif mod3p(i)<=T3(2) 
mod3(i)=l; 

elseif mod3p(i)<=T3(3) 
mod3(i)=2; 

elseif mod3p(i)<=T3(4) 

mod3(i)=3; 

else 

mod3(i)=4; 

end 

end 

figure(l) 
whitebg 
subplot(611) 
plot(theta,mod 1 p) 
axis([-90,90,-4,4]); 

title(['Foldingwave and RNS output code (sampling 

distance:ds-,num2str(ds*3),’degree)']),', num2str(l/ds),'samples in 3 degree)']) 


64 



grid 

subplot(612) 
plot(theta,mod2p) 
axis([-90,90,-4,4]); 
ylabel('Folding waves') 
grid 

subplot(613) 

plot(theta,mod3p) 

axis([-90,90,-4,4]); 

grid 

subplot(614) 
stairs(theta,modl) 
axis([-90,90,-0.2,2.2]); 
grid 

subplot(615) 
stairs(theta,mod2) 
axis([-90,90,-0.2,3-2]); 
ylabel('RNS Output Code') 
grid 

subplot(616) 

stairs(theta,mod3) 

axis([-90,90,-0.2,4.2]); 

xlabel('Input Direction of Arrival (degrees)') 
grid 

% Find the angle from the 3-bit RNS Modulus code 
mslogl 

% Change 3-bit RNS modulus code to six-bit binary code 
bin_six 

% Find the number of glitch 
glitch 

% Plots the Resolved DOA and Quantization Error as a Function of the Input Direction of 

Arrival 

figure(2) 

whitebg 

plot(theta,th,'b') 

xlabel('Input Direction of Arrival (degrees)') 
ylabel('Resolved DOA and Quantization Error (degrees)') 
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title(['Encoding Error in Quatization (glitch= ',num2str(g_count),',sampling 

distance:ds=',num2str(ds*3),' degree)'])%, ',num2str(l/ds) 5 'samples in 3 degree)']) 

axis([-90 90 -90 90]) 

grid 

print 

figure(3) 
whitebg 
subplot(211) 
plot(theta,mod 1 p) 
axis([-90,90,-4,4]); 

xlabel('Input Direction of Arrival (degrees)') 

title(['Mode 3 Foldingwave and Output code (sampling distance:ds-,num2str(ds*3),' 

degree)'])%, ',num2str(l/ds),'samples in 3 degree)']) 

grid 

subplot(212) 

stairs(theta,modl) 

axis([-90,90,-0.2,2.2]); 

xlabel('Input Direction of Arrival ( degrees)') 

grid 

save binl 

% interpolate the Encoding error 
interpl 
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%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%°/o%% 
% This program 'Treshold' calculates the original threshold levels for setting the RNS % 
% folding wave circuit in the 0.3 degree sampling distance % 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


%Program 2 
% thres.m 

global ml 
global m2 
global m3 
global fo 
global ff 
global ds 
global M 
global L 

% Calculates the Threshold Values for Modulus ml 
i=-l:2:l; 

Tl=(pi/3)*i; 

Tl=sprintf('%6.3f,Tl); 

Tl=sscanf(Tl,'%f); 

% Calculates the Threshold Values for Modulus m2 
i=-l:l; 

T2=(pi/2)*i; 

T2=sprintf('%6.3f,T2); 

T2=sscanf(T2,'%f); 

% Calculates the Threshold Values for Modulus m3 
i=-3:2:3; 

T3=(pi/5)*i; 

T3=sprintf('%6.3f,T3); 

T3=sscanf(T3,'%f); 
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%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%This program 'RNS_ANGLE' is to determine the resolved Direction Of Arrival from % 
% the RNS output in the moduli set (3,4,5) % 

%%%%%%%%%%%%%%%%%%%%%%%%%%Vo%%VoVoVo%Vo%%%%% 0 /o°/o%% 

% program 3 
% set 3-Bit RNS Code 
res=[modl' mod2' mod3']; 

% Number of Folds Over the Entire field of View for Modulus ml 
nf=fr/fo*M/ml; 

%Total Number of Quantization Levels Over the Entire Field of View 
nql=nf*ml; 

%Generates the Quantization Level Transistion Angles 
p=0:2/nql:l; 

p=[ -rot90(p(2: size(p,2)))' p]; 
trang=l/rad*asin(p); 

%Calculates the Midvalue of Each Quantization Level 
for n=l :size(p,2)-l; 

midval(n)=trang(n)+(trang(n+1 )-trang(n))/2; 
qlw(n)=midval(n)-trang(n); 
end 

if rem(size(midval,2),2)=0; 
a=midval; 
b=size(a,2)/2; 
end 

if rem(size(midval,2),2)=l; 
a=midval; 
b=(size(a,2)+l)/2; 
end 

if size(a,2)“60; 
a=[zeros(l,b) a zeros(l,b)]; 
end 

b=size(a,2)/2; 

%Calculates the Resolved Direction of Arrival 
for u=l:M/ds; 
if res(u,:)=([0,0,0]); 
th(u)=a(l); 

elseif res(u,:)=([ 1,1,1]); 
th(u)=a(2); 
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elseif res(u,:)=([2,2,2]); 
th(u)=a(3); 

elseif res(u,:)=([0,3,3]); 
th(u)=a(4); 

elseif res(u,:)==([l,0,4]); 
th(u)=a(5); 

elseif res(u,:)=([2.1,0]); 
th(u)=a(6); 

elseif res(u,:)=([0,2,l]); 
th(u)=a(7); 

elseif res(u,:)=([l,3,2]); 
th(u)=a(8); 

elseif res(u,:)=([2,0,3]); 
th(u)=a(9); 

elseif res(u,:)=([0,l,4]); 
th(u)=a(10); 

elseif res(u,:)=([J,2,0]); 
th(u)=a(ll); 

elseif res(u,:)=([2,3,l]); 
th(u)=a(12); 

elseif res(u,:)=([0,0,2]); 
th(u)=a(13); 

elseif res(u,:)=([l,l,3]); 
th(u)=a(14); 

elseif res(u,:)=([2,2,4]); 
th(u)=a(15); 

elseif res(u,:)=([0,3,0]); 
th(u)=a(16); 

elseif res(u,:)=([l,0,l]); 
th(u)=a(17); 

elseif res(u,:)=([2,1,2]); 
th(u)=a(18); 

elseif res(u,:)—([0,2,3]); 
th(u)=a(19); 

elseif res(u,:)=([l, 3,4]); 
th(u)=a(20); 

elseif res(u,:)=([2,0,0]); 
th(u)=a(21); 

elseif res(u,:)=([0,l,l]); 
th(u)=a(22); 

elseif res(u,:)=([l,2,2]); 
th(u)=a(23); 



elseif res(u,:)=([2,3,3]); 
th(u)=a(24); 

elseif res(u,:)=([0,0,4]); 
th(u)=a(25); 

elseif res(u,:)=([l,l,0]); 
th(u)=a(26); 

elseif res(u,:)=([2,2,l]); 
th(u)=a(27); 

elseif res(u,:)=([0,3,2]); 
th(u)=a(28); 

elseif res(u,:)=([l,0,3]); 
th(u)=a(29); 

elseif res(u,:)=([2,1,4]); 
th(u)=a(30); 

elseif res(u,:)=([0,2,0]); 
th(u)=a(31); 

elseif res(u, :)=([ 1,3,1]); 
th(u)=a(32); 

elseif res(u,:)=([2,0,2]); 
th(u)=a(33); 

elseif res(u,:)=([0,l,3]); 
th(u)=a(34); 

elseif res(u,:)=([l,2,4]); 
th(u)=a(35); 

elseif res(u,:)=([2,3,0]); 
th(u)=a(36); 

elseif res(u,:)=([0,0,l]); 
th(u)=a(37); 

elseif res(u,:)=([l, 1,2]); 
th(u)=a(38); 

elseif res(u,:)=([2,2,3]); 
th(u)=a(39); 

elseif res(u,:)=([0,3,4]); 
th(u)=a(40); 

elseif res(u,:)=([l,0,0]); 
th(u)=a(41); 

elseif res(u, :)=([2,1,1]); 
th(u)=a(42); 

elseif res(u,:)=([0,2,2]); 
th(u)=a(43); 

elseif res(u,:)=([l,3,3]); 
th(u)=a(44); 



elseif res(u,:)=([2,0,4]); 
th(u)=a(45); 

elseif res(u,:)—([0,1,0]); 
th(u)=a(46); 

elseif res(u,:)=([l,2,l]); 
th(u)=a(47); 

elseif res(u,:)=([2,3,2]); 
th(u)=a(48); 

elseif res(u,:)=([0,0,3]); 
th(u)=a(49); 

elseif res(u,:)=([l, 1,4]); 
th(u)=a(50); 

elseif res(u,:)=([2,2,0]); 
th(u)=a(51); 

elseif res(u,:)=([0,3,l]); 
th(u)=a(52); 

elseif res(u,:)=([l,0,2]); 
th(u)=a(53); 

elseif res(u,:)=([2,l,3]); 
th(u)=a(54); 

elseif res(u,:)=([0,2,4]); 
th(u)=a(55); 

elseif res(u,:)=([l,3,0]); 
th(u)=a(56); 

elseif res(u,:)=([2,0,l]); 
th(u)=a(57); 

elseif res(u,:)=([0,l,2]); 
th(u)=a(58); 

elseif res(u,:)=([l,2,3]); 
th(u)=a(59); 
else res(u, :)=([2,3,4]) ; 
th(u)=a(60); 
end 
end 
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%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%°/o%% 
% This program 'SIX_BINARY' converts the RNS code output % 

% to 6-bit binary code at the instant input angle % 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 0 /o 0 /o 

% program 4 
% bin_six.m 
% make six-binary code 
for u=l :M/ds 
if res(u,:)=([0,0,0]) 
thl(u)=l; 

binary(u,:)=([0 0 0 0 0 0]); 

elseif res(u,:)=([l,l,l]) 
thl(u)=2; 

binary(u,:)=([0 0 0 0 0 1]); 

elseif res(u,:)=([2,2,2]) 
thl(u)=3; 

binaiy(u,:)=([0 0 0 01 0]); 

elseif res(u,:)=([0,3,3]) 
thl(u)=4; 

binary(u,:)=([0 0 0 0 1 1]); 

elseif res(u,:)==([l > 0,4]) 
thl(u)=5; 

binary(u,:)=([0 0 0 1 0 0]); 

elseif res(u,:)==([2,1,0]) 
thl(u)=6; 

binary(u,:)=([0 0 0 1 0 1]); 

elseif res(u,:)=( [0,2,1]) 
thl(u)=7; 

binary(u,:)=([0 0 0 1 1 0]); 

elseif res(u,:)=([l ,3,2]) 
thl(u)=8; 

binary(u,:)=([0 0 0 1 1 1]); 
elseif res(u,:)=([2,0,3]) 
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thl(u)=9; 

binary(u,:)=([0 01 0 0 0]); 

elseif res(u,:)=([0,l,4]) 
thl(u)=10; 

binary(u,:)=([0 01 0 01]); 

elseif res(u,:)=([l ,2,0]) 
thl(u)=ll; 

binary(u,:)=([0 0101 0]); 

elseif res(u,:)=([2,3,l]) 
thl(u)=12; 

binary(u,:)=([0 0101 1]); 

elseif res(u,:)=([0,0,2]) 
thl(u)=13; 

binary(u,:)=([0 01 1 0 0]); 

elseif res(u, :)=([ 1,1,3]) 
thl(u)=14; 

binary(u,:)=([0 0 110 1]); 

elseif res(u,:)=([2,2,4]) 
thl(u)=15; 

binary(u,:)=([0 0 111 0]); 

elseif res(u,:)=([0,3,0]) 
thl(u)=16; 

binary(u,:)=([0 0 111 1]); 

elseif res(u,:)=([l,0,l]) 
thl(u)=17; 

binary(u,:)=([0 1 0 0 0 0]); 

elseif res(u,:)=([2,l,2]) 
thl(u)=18; 

binary(u,:)=([0 1 0 0 01]); 

elseif res(u,:)=([0,2,3]) 
thl(u)=19; 

binary(u,:)=([0 10 01 0]); 



elseif res(u, :)==([ 1,3,4]) 
thl(u)=20; 

binary(u,:)=([0 10 01 1]); 

elseif res(u,:)=([2,0,0]) 
thl(u)=21; 

binary(u,:)=([0 1010 0]); 

elseif res(u, :)=([0,1,1]) 
thl(u)=22; 

binary(u,:)=([0 1010 1]); 

elseif res(u,:)=([l,2,2]) 
thl(u)=23; 

binaiy(u,:)=([0 1011 0]); 

elseif res(u,:)={[2,3,3]) 
thl(u)=24; 

binary(u,:)=([0 101 1 1]); 

elseif res(u,:)=([0,0,4]) 
thl(u)=25; 

binary (u,:)=([0 1 1 0 0 0]); 

elseif res(u,:)=([l, 1,0]) 
thl(u)=26; 

binary(u,:)=([0 1 10 01]); 

elseif res(u,:)=([2,2,l]) 
thl(u)=27; 

binary(u,:)=([0 1101 0]); 

elseif res(u,:)=([0,3,2]) 
thl(u)=28; 

binary(u,:)=([0 110 1 1]); 

elseif res(n,:)=([l,0,3]) 
thl(u)=29; 

binary(u,:)=([0 1110 0]); 


elseif res(u,:)=([2,l,4]) 



thl(u)=30; 

binary(u,:)=([0 1 1 101]); 

elseif res(u,:)=([0,2,0]) 
thl(u)=31; 

binary(u,:)=([0 1111 0]); 

elseif res(u,:)=([ 1,3,1]) 
thl(u)=32; 

binary(u,:)=([0 1111 1]); 

elseif res(u,:)=([2,0,2]) 
thl(u)=33; 

binary(u,:)=([l 0 0 0 0 0]); 

elseif res(u,:)=([0,l,3]) 
thl(u)=34; 

binary(u,:)=([l 0 0 0 0 1]); 

elseif res(u,:)=([l,2,4]) 
thl(u)=35; 

binary(u,:)=([l 0 0 0 1 0]); 

elseif res(u,:)=([2,3,0]) 
thl(u)=36; 

binaiy(u,:)=([l 0 0 0 1 1]); 

elseif res(u,:)=([0,0,l]) 
thl(u)=37; 

binary(u,:)=([l 0 01 0 0]); 

elseif res(u, :)=([1,1,2]) 
thl(u)=38; 

binary(u,:)=([l 0 010 1]); 

elseif res(u,:)=([2,2,3]) 
thl(u)=39; 

binary(u,:)=([l 0 011 0]); 

elseif res(u,:)=([0,3,4]) 
thl(u)=40; 

binary(u,:)=([! 0 011 1]); 



elseif res(u,:)=([ 1,0,0]) 
thl(u)=41; 

binary(u,:)=([l 0 1 0 0 0]); 

elseif res(u, :)=([2,1,1]) 
thl(u)=42; 

binary(u,:)=([l 010 0 1]); 

elseif res(u,:)=([0,2,2]) 
thl(u)=43; 

binary(u,:)=([l 0101 0]); 

elseif res(u,:)=([l,3,3]) 
thl(u)=44; 

binaiy(u,:)=([l 0101 1]); 

elseif res(u,:)=([2,0,4]) 
thl(u)=45; 

binary(u,:)=([l 0110 0]); 

elseif res(u,:)=([0,l,0]) 
thl(u)=46; 

binary(u,:)=([l 01 101]); 

elseif res(u,:)=([l,2,l]) 
thl(u)=47; 

binary(u,:)=([l 0 111 0]); 

elseif res(u,:)=([2,3,2]) 
thl(u)=48; 

binary(u,:)=([l 0 1111]); 

elseif res(u,:)=([0,0,3]) 
thl(u)=49; 

binary(u,:)=([l 1 0 0 0 0]); 

elseif res(u, :)==([ 1,1,4]) 
thl(u)=50; 

binary(u,:)=([! 1 0 0 0 1]); 


elseif res(u,:)=([2,2,0]) 



thl(u)=51; 

binary(u,:)=([l 10 01 0]); 

elseif res(u,:)=([0,3,l]) 
thl(u)=52; 

binary(u,:)=([l 10 01 1]); 

elseif res(u,:)=([l ,0,2]) 
thl(u)=53; 

binary(u,:)=([l 1010 0]); 

elseif res(u,:)=([2,l,3]) 
thl(u)=54; 

binary(u,:)=([! 1010 1]); 

elseif res(u,:)=([0,2,4]) 
thl(u)=55; 

binary(u,:)=([l 1011 0]); 

elseif res(u,:)=([l ,3,0]) 
thl(u)=56; 

binary(u,:)=([l 10111]); 

elseif res(u,:)=([2,0,l]) 
thl(u)=57; 

binary(u,:)=([! 1 1 0 0 0]); 

elseif res(u,:)=([0,l,2]) 
thl(u)=58; 

binary(u,:)=([l 110 0 1]); 

elseif res(u,:)=([l,2,3]) 
thl(u)=59; 

binary(u,:)=([l 110 1 0]); 

elseif res(u,:)=([2,3,4]) 
thl(u)=60; 

binary(u,:)=([l 1 10 1 1]); 

end 

end 

%binary' 



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%°/o 
% This program 'Shift_LSB' performs the interpolation using Shift LSB Method % 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% program 5 
% interpl.m 

% set the six-addittional threshol 
adth3 

% This is the operating loop until 'M/ds' number of parity samples 
i=l; 

while i<M/ds 

a=l; % initial value of'a' is 1 

% number of continuous error parity in 5 parity samples 
% for example, when error is 1, a=2 
while parity (i)~=l % to find the first T parity(error parity) 
i=i+l; % i is the location of error 

end 

while parity(i)~=0 % to find the next error parity('O'value) 
i=i+l; 
end 

while parity (i)~=l % to find the continuous number of error parity 
i=i+l; 
a=a+l; 

if i>=M/ds+1,break,end % if the number of parity is 
% bigger than M/ds,then finish loop ! 
end 

if i>=M/ds+1,break,end 
i=i-a+l; % real location value of'O'parity 
b=l; % highest level of LSB value 
while binary(i-1 ,b)=binary(i+a-1 ,b) 
b=b+l; 
if b=7 
break 
end 
end 
if b=l 
for k=i+3:-l:i 
forj=l:6 

binary(k,j )=binary(k-1 ,j); 
end 
end 
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end 

% when the number of parity error is one (a=2) 
if a=2 

% change the MSB value 
forj=l:b-l 

binary(i+3j)=binary(i+2,j); 
binary(i+2 j)=binary(i+l j); 
binary (i+1 ,j )=binary(i+1 j); 
binary( i, j)=binary(i-l j); 
end 

% change the LSB value 
if b<7 

fork=i+3:-l:i 

forj=b:6 

binary(k j)=binary(k-l ,j); 
end 
end 
end 

% when the number of parity error is two (a=3) 
elseif a=3 

% change the MSB value 
forj=l:b-l 

binary(i+3 ,j)=binary(i+2,j); 
binary(i+2j)=binary(i-l j); 
binary(i+1 ,j)=binary(i-1 j); 
binary( i,j)=binary(i-l,j); 

end 

% change the LSB value 
if b<7 

for k=i+3:-l:i 
forj=b:6 

binary(k j)=binary(k-l ,j); 
end 
end 
end 

% when the number of parity error is three (a=4) 
elseif a=4 

% change the MSB value 
for j=l:b-l 
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binary (i-t-3 j )=binary (i-1 j); 
binary(i+2 j)=binaiy(i-1 j); 
binary(i+l ,j)=binary(i-l j); 
binary( i, j)=binary(i-lj); 
end 

% change the LSB value 
ifb<7 

for k=i+3:-l:i 
forj=b:6 

binary (k,j )=binary (k-1 j); 
end 
end 
end 

% when the number of parity error is four 
%a=5 

% change the all binary value 
else 

for k=i+3:-l :i+l 
for j=l:6 

binary (k,j )=binary(k-1 j); 

end 

end 

end % if a=2 

% printout 6-bit binary-code after interpolation 
if (i+3) < (M/ds) 
v=i-l:i+3; 
parity(v) 
binary (v,:)' 
end 
i=i+a; 

end % while i<=M/ds 

% change 6-bit binary code to Modulus (ml,m2,m3) 
antbin 

% Calculates the Resolved Direction of Arrival after interpolation 
mslog2 

% find the number of glitch 

gi_count=0; 

for u=l :M/ds-l 
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max(u+1 )=th2(u+1 )-th2(u); 
if u<M/(2*ds) 
ifmax(u+l)> 11 

parity(u+l )=0; % possible bad samples 

gi_count=gi_count+l; 
u+l,max(u+l) 
parity(u:u+4) 
binary(u:u+4,:)' 
else 

parity(u+l)=l; % sample is good 
end 
else 

if max(u+l)> 11 

parity(u)=0; % possible bad samples 
gi_count=gi_count+l; 
u+l,max(u+l) 
parity(u-l :u+3) 
binary(u-l :u+3,:)' 
else 

# 

parity(u)=l; % sample is good 
end 

end% ifu<M/(2*ds) 
end % for 


figure(5) 

whitebg 

plot(theta,th2,'b') 

xlabel('Input Direction of Arrival (degrees)') 

ylabel('Resolved DOA and Quantization Error ( degrees)') 

title(['Interpolation Result,',num2str(LSB_percent),' % LSB-Shift Method ( glitch= 

',num2str(gi_count),', sampling distance:ds- ,num2str(ds*3),' degree)']) 

axis([-90 90 -90 90]) 

save bin2 Tl_add gi_count binary parity; 
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%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
% This progr am 'ADD_THRESH' calculates the additional threshold levels % 

% to remove the quantization error % 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 0 /o%%%%% 0 /o 0 /o 0 /o 0 /o 0 /o°/o 


% program 6 
% adth3.m 

global ml 
global m2 
global m3 
global fo 
global fr 
global ds 
global rad 
global M 
global theta 
global thet 

%LSBjpercent =input('What is the percentage of LSB ?(%)') 
incre=0.6*LSB_percent/l00; % Set the threshold by 1 degree 

% set the additional thresholds 
T1 _add( 1 )=2*pi* incre/9-pi; 

Tl_add(2)=2*pi*(3-incre)/9-pi; 

T1 _add(3 )=2 *pi * (3+incre)/9-pi; 

T1 _add(4)=2 *pi* (6-incre)/9-pi; 
Tl_add(5)=2*pi*(6+incre)/9-pi; 
Tl_add(6)=2*pi*(9-incre)/9-pi; 

Tl_add 

% calculate the six-comparator output 
for i=l:M/ds 

if modlp(i) <= Tl_add(l) 
count(i) = 0; % Bad sample 

elseif modlp(i) <= Tl_add(2) 
count(i) =1; % Good sample 

elseif modlp(i) <= Tl_add(3) 
count(i) = 2; 

elseif modlp(i) <= Tl_add(4) 
count(i) - 3; 

elseif modlp(i) <= Tl_add(5) 
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count(i) = 4; 

elseif modlp(i) <= Tl_add(6) 
count(i) = 5; 
else 

count(i) = 6; 
end 
end 


% parity cicuit 
for i=l:M/ds 
if rem(count(i),2) — 0 
parity(i)=0; % possible bad 
else 

parity(i)=1; % good parity 

end 
end 

figure(4) 
whitebg 
subplot(211) 
plot(theta,modlp) 
axis([-90,90,-4,4]); 

xlabel('Input Direction of Arrival ( degrees)') 
grid 

subplot(212) 
stairs(thet,count) 
axis([-90,90,-0.2,6.2]); 

xlabel('Input Direction of Arrival ( degrees)') 

grid 

print 
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%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%°/o%%% 
% This program 'ANT_BINARY' converts the six-bit binary code % 

% to the RNS code after interpolation the encoding error % 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% program 7 
% antbin.m 

% change six-binary code to RNS Code 
foru=l:M/ds 

if binary(u,:)=([0 0 0 0 0 0]) 
res2(u,:)=([0 ; 0,0]); 
thl(u)=l; 

elseif binary(u, :)=([0,0,0,0,0,1 ]) 

res2(u,:)=([l,l,l]); 

thl(u)=2; 

%binary(u,:) 

%res(u,:) 

elseif binary(u,:)—([0 0 0 0 1 0]) 
res2(u,:)=([2 5 2,2]); 
thl(u)=3; 

elseif binary(u,: )=( [0 0 0 01 1]) 
res2(u,:)=([0 5 3,3]); 
thl(u)=4; 

elseif binary(u,:)=([0 0 01 0 0]) 
res2(u,:)=([l ,0,4]); 
thl(u)=5; 

elseif binaiy(u,:) : =([0 0 01 0 1]) 
res2(u,:)=([2,l,0]); 
thl(u)=6; 

elseif binary(u,:)=([0 0 01 1 0]) 
res2(u,:)=([0,2,l]); 
thl(u)=7; 

elseif binary(u,:)==([0 0 01 1 1]); 
res2(u,:)=([l ,3,2]); 
thl(u)=8; 

elseif binary(u,:)=([0 01 0 0 0]); 

res2(u,:)=([2,0,3]); 

thl(u)=9; 

elseif binary(u,:)=([0 0 1 0 0 1]); 
res2(u,:)=([0,l,4]); 
thl(u)=10; 

elseif binary(u,:)=([0 0 10 1 0]); 
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res2(u,:)=([l,2,0]); 

thl(u)=ll; 

elseif binary(u,:)==([0 0101 1]); 

res2(u,:)=([2,3,l]); 

thl(u)=12; 

elseif binary(u,:)=([0 01 10 0]); 
res2(u,:)=([0,0,2]); 
thl(u)=13; 

elseif binary(u,:)==([0 01 101]); 
res2(u,:)=([l,l,3]); 
thl(u)=14; 

elseif binary(u,:)=([0 0111 0]); 
res2(u,:)=([2,2,4]); 
thl(u)=15; 

elseif binary(u,:)=([0 01111]); 
res2(u,:)=([0,3,0]); 
thl(u)=16; 

elseif binary(u,:)=([0 1 0 0 0 0]); 
res2(u,:)=([ 1,0,1]); 
thl(u)=17; 

elseif binary(u,:)=([0 1 0 0 0 1]); 
res2(u,:)=([2,l,2]); 
thl(u)=18; 

elseif binary(u,:)==([0 10 01 0]); 
res2(u,:)=([0,2,3]); 
thl(u)=19; 

elseif binary(u,:)—([0 10 01 1]); 
res2(u,:)=([l,3,4]); 
thl(u)=20; 

elseif binary(u,:)==([0 10 10 0]); 
res2(u,:)=([2,0,0]); 
thl(u)=21; 

elseif binary(u,:)=([0 10 10 1]); 
res2(u,:)=([0,l,l]); 
thl(u)=22; 

elseif binary(u,:)=([0 1011 0]); 
res2(u,:)=([l,2,2]); 
thl(u)=23; 

elseif binary^O^CtO 10 11 1]); 
res2(u,:)=([2,3,3]); 
thl(u)=24; 

elseif binary(u,:)=([0 1 1 0 0 0]); 



res2(u,:)=([0,0,4]); 

thl(u)=25; 

elseif binary(u,:)==([0 110 0 1]); 
res2(u,:)=([U,0]); 
thl(u)=26; 

elseif binary(u,:)=([0 1101 0]); 
res2(u,:)=([2,2,l]); 
thl(u)=27; 

elseif binaiy(u,:)=([0 1101 1]); 
res2(u,:)=([0,3,2]); 
thl(u)=28; 

elseif binary(u,:)==([0 1110 0]); 
res2(u,:)=([l,0,3]); 
thl(u)=29; 

elseif binary(u,:)=([0 1110 1]); 
res2(u,:)=([2,l,4]); 
thl(u)=30; 

elseif binary(u,:)==([0 1111 0]); 
res2(u,:)=([0,2,0]); 
thl(u)=31; 

elseif binary(u,:)=([0 1111 1]); 

res2(u,:)=([l,3,l]); 

thl(u)=32; 

elseif binary(u,:)=([l 0 0 0 0 0]); 
res2(u,:)=([2,0,2]); 
thl(u)=33; 

elseif binary(u,:)=([l 0 0 0 0 1]); 
res2(u,:)=([0,l,3]); 
thl(u)=34; 

elseif binary(u,:)=([l 0 0 0 1 0]); 
res2(u,:)=([l,2,4]); 
thl(u)=35; 

elseif binary(u,:)=([l 0 0 0 1 1]); 
res2(u,:)=([2,3,0]); 
thl(u)=36; 

elseif binary(u,:)==([l 0 0 1 0 0]); 
res2(u,:)=([0,0,l]); 
thl(u)=37; 

elseif binary(u,:)=([l 0 010 1]); 
res2(u,:)=([l,l,2]); 
thl(u)=38; 

elseif binary(u,:)=([! 0 0 11 0]); 



res2(u,:)=([2,2,3]); 

thl(u)=39; 

elseif binary(u,:)==([l 0 011 1]); 
res2(u,:)=([0,3,4]); 
thl(u)=40; 

elseif binary(u,:)==([l 01 0 0 0]); 
res2(u,:)=([! AO]); 
thl(u)=41; 

elseif binary(u,:)==([l 010 0 1]); 
res2(u,:)=([2,l,l]); 
thl(u)=42; 

elseif binary(u, :)==([ 1 0101 0]); 

res2(u,:)=([0,2,2]); 

thl(u)=43; 

elseif binary(u, :)==([ 1 0101 1]); 

res2(u,:)=([13,3]); 

thl(u)=44; 

elseif binary(u,:)=([l 0110 0]); 
res2(u,:)=([2,0,4]); 
thl(u)=45; 

elseif bmary(u,:)==([l 01 101]); 
res2(u,:)=([0,l,0]); 
thl(u)=46; 

elseif binary(u,:)=([l 0111 0]); 
res2(u,:)=([l,2,l]); 
thl(u)=47; 

elseif binaiy(u,:)=([l 0111 1]); 
res2(u,:)=([2,3,2]); 
thl(u)=48; 

elseif binary(u,:)=([l 1 0 0 0 0]); 
res2(u,:)=([0,0,3]); 
thl(u)=49; 

elseif binaiy(u,:)=([l 10 0 0 1]); 
res2(u, :)=([1,1,4]); 
thl(u)=50; 

elseif binary(u,:)=([l 10 01 0]); 
res2(u,:H[2,2,0]); 
thl(u)=51; 

elseif binary(u,:)=([l 10 01 1]); 
res2(u, :>=([0,3,1 ]); 
thl(u)=52; 

elseif binary(u,:)=([! 10 10 0]); 



res2(u, :)=([!,0,2]); 
thl(u)=53; 

elseif binary (u,:)=([l 10 10 1]); 
res2(u,:)=([2,l,3]); 
thl(u)=54; 

elseif binary(u,:)=([l 1011 0]); 
res2(u,:)=([0,2,4]); 
thl(u)=55; 

elseif binary(u, :)==([ 1 1011 1]); 
res2(u,:)=([l,3,0]); 
thl(u)=56; 

elseif binary(u,:)=([l 1 1 0 0 0]); 
res2(u,:)=([2,0,l]); 
thl(u)=57; 

elseif binary(u,:)=([l 110 0 1]); 
res2(u,:)=([0,l,2]); 
thl(u)=58; 

elseif binary(u,:)=([l 110 1 0]); 
res2(u,:)=([l ,2,3]); 
thl(u)=59; 

elseif binary(u,:)=([l 110 1 1]); 
res2(u,:)=([2,3,4]); 
thl(u)=60; 
end 
end 


%binary' 



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 0 /o°/o% 

% This program 'RNS_angle2' is to determine the resolved Direction Of Arrival % 
% the RNS output in the moduli set (3,4,5) % 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 0 / 0 % 0 / 0%%%%% 0 /0 


% Program 8 
% mslog2.m 

% Number of Folds Over the Entire field of View for Modulus ml 
nf=fr/fo*M/ml; 

%Total Number of Quantization Levels Over the Entire Field of View 
nql=nf*ml; 

%Generates the Quantization Level Transition Angles 
p=0:2/nql:l; 

p=[ -rot90(p(2:size(p,2)))' p]; 
trang=l/rad*asin(p); 

%Calculates the Midvalue of Each Quantization Level 
forn=l:size(p,2)-l 

midval(n)=trang(n)+(trang(n+l)-trang(n))/2; 

%qlw(n)=midval(n)-trang(n); 

end 

if rem(size(midval,2),2)=0 

a=midval; 

b=size(a,2)/2; 

end 

if rem(size(midval,2),2)=l 

a=midval; 

b=(size(a,2)+l)/2; 

end 

if size(a,2)~=60 
a=[zeros(l,b) a zeros(l,b)]; 
end 

b=size(a,2)/2; 

%Calculates the Resolved Direction of Arrival 
for u=l:M/ds 
if res2(u,:)=([0,0,0]); 
th2(u)=a(l); 

elseif res2(u,:)=([l ,1,1]); 
th2(u)=a(2); 

elseif res2(u,:)=([2,2,2]); 
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th2(u)=a(3); 

elseif res2(u,:)=([0,3,3]); 
th2(u)=a(4); 

elseif res2(u,:)=([l,0,4]); 
th2(u)=a(5); 

elseif res2(u,:)—([2,1,0]); 
th2(u)=a(6); 

elseif res2(u, :)=([0,2,1 ]); 
th2(u)=a(7); 

elseif res2(u, :)=([ 1,3,2]); 
th2(u)=a(8); 

elseif res2(u,:)=([2,0,3]); 
th2(u)=a(9); 

elseif res2(u,:)=([0,l,4]); 

th2(u)=a(10); 

elseif res2(u,:)=([ 1,2,0]); 

th2(u)=a(l 1); 

elseif res2(u,:)=([2,3,l]); 

th2(u)=a(12); 

elseif res2(u,:)=([0,0,2]); 

th2(u)=a(13); 

elseif res2(u,:)=([l,l,3]); 

th2(u)=a(14); 

elseif res2(u,:)=([2,2,4]); 

th2(u)=a(15); 

elseif res2(u,:)=([0,3,0]); 

th2(u)=a(16); 

elseif res2(u,:)=([l,0,l]); 

th2(u)=a(17); 

elseif res2(u,:)=([2,1,2]); 

th2(u)=a(18); 

elseif res2(u,:)=([0,2,3]); 

th2(u)=a(19); 

elseif res2(u,:)=([l,3,4]); 

th2(u)=a(20); 

elseif res2(u,:)=([2,0,0]); 

th2(u)=a(21); 

elseif res2(u,:)=([0,l,l]); 

th2(u)=a(22); 

elseif res2(u,:)=([ 1,2,2]); 

th2(u)=a(23); 

elseif res2(u,:)=([2,3,3]); 



th2(u)=a(24); 

elseif res2(u,:)=([0,0,4]); 

th2(u)=a(25); 

elseif res2(u,:)=([l,l,0]); 

th2(u)=a(26); 

elseif res2(u,:)=([2,2,l]); 

th2(u)=a(27); 

elseif res2(u,:)=([0,3,2]); 

th2(u)=a(28); 

elseif res2(u,:)=([ 1,0,3]); 

th2(u)=a(29); 

elseif res2(u,:)=([2,l,4]); 

th2(u)=a(30); 

elseif res2(u,:)—([0,2,0]); 

th2(u)=a(31); 

elseif res2(u,:)=([l,3,l]); 

th2(u)=a(32); 

elseif res2(u,:)=([2,0,2]); 

th2(u)=a(33); 

elseif res2(u,:)=([0,l,3]); 

th2(u)=a(34); 

elseif res2(u,:)—([1,2,4]); 

th2(u)=a(35); 

elseif res2(u,:)=([2,3,0]); 

th2(u)=a(36); 

elseif res2(u,:)=([0,0,l]); 

th2(u)=a(37); 

elseif res2(u,:)==([l, 1,2]); 

th2(u)=a(38); 

elseif res2(u,:)=([2,2,3]); 

th2(u)=a(39); 

elseif res2(u,:)=([0,3,4]); 

th2(u)=a(40); 

elseif res2(u,:)=([l,0,0]); 

th2(u)=a(41); 

elseif res2(u,:)=([2,1,1 ]); 

th2(u)=a(42); 

elseif res2(u,:)=([0,2,2]); 

th2(u)=a(43); 

elseif res2(u,:)=([l,3,3]); 

th2(u)=a(44); 

elseif res2(u,:)=([2,0,4]); 



th2(u)=a(45); 
elseif res2(u,:)=([0,l,0]); 
th2(u)=a(46) ; 
elseif res2(u, :)=([ 1,2,1 ]); 
th2(u)=a(47); 
elseif res2(u,:)=([2,3,2]); 
th2(u)=a(48); 
elseif res2(u,:)=([0,0,3]); 
th2(u)=a(49); 
elseif res2(u,:)=([l,1,4]); 
th2(u)=a(50); 
elseif res2(u,:)=([2,2,0]); 
th2(u)=a(51); 
elseif res2(u,:)=([0,3,l]); 
th2(u)=a(52); 
elseif res2(u,:)=([l,0,2]); 
th2(u)=a(53); 
elseif res2(u,:)=([2,l ,3]); 
th2(u)=a(54); 
elseif res2(u,:)=([0,2,4]); 
th2(u)=a(55); 
elseif res2(u,:)=([ 1,3,0]); 
th2(u)=a(56); 
elseif res2(u,:)==([2,0,l]); 
th2(u)=a(57); 
elseif res2(u,:)=([0,1,2]); 
th2(u)=a(58); 
elseif res2(u,:)=([l,2,3]); 
th2(u)=a(59); 
else res2(u,:)=([2,3,4]); 
th2(u)=a(60); 
end 
end 



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 0 /o% 0 /o%%%% 
% This program 'Last_Good' performs the interpolation using Last good-sample shift % 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 0 /o 0 /o 0 /o 0 /o%%% 

% program 9 
% int_copy 

% set the six-addittional threshol 
adth3 

% This is the operating loop until 'M/ds' number of parity samples 

i=l; 

while i<M/ds 

a=l; % initial value of'a' is 1 

% number of continuous error parity in 5 parity samples 
% for example, when error is 1, a=2 
while parity(i)~=l % to find the first T parity(error parity) 
i=i+l; % i is the location of error 

end 

while parity(i)~=0 % to find the next error parity('O'value) 
i=i+l; 
end 

while parity(i)~=l % to find the continuous number of error parity 
i=i+l; 
a=a+l; 

if i>=M/ds+l,break,end % if the number of parity is 
% bigger than M/ds,then finish loop ! 
end 

if i>=M/ds+l,break,end 

i=i-a+l; % real location value of'O'parity 

b= 1; % highest level of LSB value 

% Copy the last good sample to error samples 
for k=i+a:-l :i 
for j=l:6 

binary (k,j)=binary(i-2j); 
end 
end 
i=i+l; 

end % while i<=M/ds 

% change 6-bit binary code to Modulus (ml,m2,m3) 
ant bin 
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% Calculates the Resolved Direction of Arrival after interpolation 
mslog2 

% find the number of glitch 
%parity=ones( 1 ,M/ds); 
gi_count=0; 
for u=l :M/ds-l 
max(u+1 )=th2(u+1 )-th2(u); 
if u<M/(2*ds) 
if max(u+l)> 11 

parity(u+l)=0; % possible bad samples 
gi_count=gi_count+l; 
u+l,max(u+l) 
parity(u:u+4) 
binary(u:u+4,:)' 
else 

parity(u+l)=l; % sample is good 
end 
else 

ifmax(u+l)> 11 

parity(u)=0; % possible bad samples 
gi_count=gi_count+l; 
u+l,max(u+l) 
parity(u-l :u+3) 
binary(u-l:u+3,:)' 
else 

parity(u)=l; % sample is good 
end 

end% ifu<M/(2*ds) 
end % for 

figure(5) 

whitebg 

plot(theta,th2,'b') 

xlabel('Input Direction of Arrival (degrees)') 
ylabel('Resolved DOA and Quantization Error (degrees )') 

title(['Interpolation Result,',num2str(LSBjercent),' % Shift Last-Good Bit Method (glitch= 
%',num2str(gi_count),', sampling distance:ds- ,num2str(ds*3),' degree)']) 
axis([-90 90 -90 90]) 
print 

save bin3 Tl_add gi_count binary parity; 



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%°/«% %%0/ ° %%% 
% This program 'RANDJLSB' performs the interpolation using random LSB Method 

% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 0 /o 0 /o%%°/o%%%%%% %0 /o % 

%program 10 
% int_rand 

% set the six-addittional threshol 
adth3 

% This is the operating loop until 'M/ds' number of parity samples 

i=l; 

while i<M/ds 

a=l; % initial value of'a' is 1 

% number of continuous error parity in 5 parity samples 
% for example, when error is 1, a=2 

while parity(i)~=l % to find the first T parity(error parity) 
i=i+l; % i is the location of error 

end 

while parity(i)~=0 % to find the next error parity('0'value) 
i=i+l; 
end 

while parity(i)~=T % to find the continuous number of error parity 
i=i+l; 
a=a+l; 

if i>=M/ds+l,break,end % if the number of parity is 
% bigger than M/ds,then finish loop ! 

end 

if i>=M/ds+1 ,break,end 
i=i-a+l; % real location value of'O'parity 
b=l; % highest level of LSB value 
while binary(i-l,b)=binary(i+a-l,b) 
b=b+l; 
ifb=7 
break 
end 
end 
if b=l 
for k=i+3:-l:i 

binary(k,:)=rem(randperm(6),2); 

end 
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end 


% when the number of parity error is one (a=2) 
ifa=2 

% change the MSB value 
forj=l:b-l 

binary(i+3j)=binary(i+2j); 
binary (i+2j )=binary(i+1 ,j); 
binary(i+l j)=binary(i+l j); 
binary( i, j)=binary(i-l j); 
end 

% randomize the LSB value 
if b<7 

fork=i+3:-l:i 

binary (k, :)=rem(randperm(6),2); 
end 
end 

% when the number of parity error is two (a=3) 
elseif a=3 

% change the MSB value 
forj=l:b-l 

binary(i+3 j)=binary(i+2 j); 
binary(i+2j)=binary(i-l,j); 
binary(i+l j)=binary(i-1 j); 
binary( i, j)=binary(i-l j); 
end 

% change the LSB value 
if b<7 

for k=i+3:-l:i 
forj=b:6 

binary(k j )=binary(k-1 j); 
end 
end 
end 

% when the number of parity error is three (a=4) 
elseif a=4 

% change the MSB value 
for j=l:b-l 

binary(i+3 ,j)=binary(i-1 j); 
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binary(i+2 j)=binary(i-l j); 
binary(i+lj)=binary(i-lj); 
binary( i, j>=binary(i-l j); 
end 

% change the LSB value 
ifb<7 

fork=i+3:-l:i 

binary(k,:)=rem(randperm(6),2); 

end 

end 

% when the number of parity error is four 
% a=5 

% change the all binary value 
else 

for k=i+3:-l:i+l 
forj=l:6 

binary(kj )=binary (k-1 j); 
end 
end 

end %ifa=2 

% printout 6-bit binary-code after interpolation 
i 

if (i+3) < (M/ds) 
v=i-l:i+3; 
parity(v) 
binary(v,:)' 
end 
i=i+a; 

end % while i<=M/ds 

% change 6-bit binary code to Modulus (ml,m2,m3) 
ant_bin 

% Calculates the Resolved Direction of Arrival after interpolation 
mslog2 

% find the number of glitch 
parity=ones( 1 ,M/ds); 
gi_count=0; 
for u=l :M/ds-l 
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max(u+1 )=th2(u+1 )-th2(u); 
if u<M/(2*ds) 
ifmax(u+l)> 11 

parity(u+l)=0; % possible bad samples 

gi_count=gi_count+l; 
u+l,max(u+l) 
parity(u:u+4) 
binary(u:u+4,:)' 
else 

parity(u+l)=l; % sample is good 
end 
else 

if max(u+l)> 11 

parity(u)=0; % possible bad samples 
gi_count=gi_count+l; 
u+l,max(u+l) 
parity(u-l :u+3) 
binary (u-l:u+3,:)' 
else 

parity(u)=l; % sample is good 

end 

end% if u<M/(2*ds) 
end % for 

figure(5) 

whitebg 

plot(theta,th2,'b') 

xlabel('Input Direction of Arrival (degrees )') 
ylabel('Resolved DOA and Quantization Error (degrees )') 

title(['Interpolation Result,',num2str(LSB_percent),' % Random-LSB Method ( glitch= 
%',num2str(gi_count),', sampling distance:ds-,num2str(ds*3),' degree)']) 
axis([-90 90 -90 90]) 

%print 

save bin4 Tl_add gi_count binary parity; 
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